The spectra of mixed 3 He- 4 He droplets 
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The diffusion Monte Carlo technique is used to calculate and analyze the excitation spectrum of 
3 He atoms bound to a cluster of 4 He atoms, by using a previously determined optimum filling of 
single-fermion orbits with well defined orbital angular momentum L, spin S and parity quantum 
numbers. The study concentrates on the energies and shapes of the three kinds of states for which 
the fermionic part of the wave function is a single Slater determinant: maximum L or maximum 
S states within a given orbit, and fully polarized clusters. The picture that emerges is that of 
systems with strong shell effects whose binding and excitation energies are essentially determined 
by averages over configuration at fixed number of particles and spin, i.e., by the monopole properties 
of an effective Hamiltonian. 

PACS numbers: 36.40.-c 61.46.+W 



I. INTRODUCTION 

The study of liquid helium in confined geometries is 
currently an active area of experimental and theoretical 
researches Helium droplets are weakly bound quantum 
systems, as a consequence of their small atomic mass and 
the particular form of the associated van der Waals in- 
teraction. For 4 Hc, clusters are bound for any number of 
atoms, while for 3 Hc, it takes about 30 — 32 atoms 3 ^ to 
form a bound system, due to the larger zero-point motion 
and the Pauli Principle. The case of mixed 3 He- 4 He clus- 
ters is very interesting, since they are made of particles 
with different statistics and masses interacting through 
the same potential. The theoretical calculations predict 
the existence of instability islands for a sufficiently small 
number of 4 He atomsASiLS 

Experimentally, small helium clusters are produced by 
free jet expansion of the gas. Their mass is then measured 
by diffraction through a transmission grating followed by 
a mass spectrometer detector^ Pure 4 He clusters, con- 
taining up to eight atoms, and mixed clusters, containing 
one 3 He and up to six 4 He atoms, have been detected us- 
ing a grating with 100 nm period^ and even the very 
weakly bound dimer 4 He2 has been unambiguously de- 
tected^ The experimental setup of Ref. |9| has been im- 
proved to detect droplet sizes up to 25-30 amu, and very 
small mixed systems have been definitely identified^ 

There have been several theoretical studies of a single 
3 He atom in a medium size 4 He cluster, either through 
a density functional approach^ 3 *^ or microscopic meth- 
oc j Si i5 4 i6 i i7 ah 4 HeAr fl 3 He clusters form bound states 
for Nb > 2. The excess in kinetic energy pushes the 



3 He atom to the surface, resulting in a a quasi two- 
dimensional wave function similar to the Andreev state 
describing one 3 He impurity in 4 He bulk. Recently, the 
ordering of the single-particle orbital states has been es- 
tablished^ Larger combinations of 3 He and 4 He atoms 
have been studied employing a non-local finite range den- 
sity functional4^2i2& Previous variational microscopic 
studies of these mixed systems have been carried out for 
the 4 He2 3 He2 cluster^ and for clusters with Nb < 8, 
Np < 20,« More recently, the diffusion Monte Carlo 
(DMC) method has been applied to droplets with Np < 
3, Nb < 17 in Refs. 00 and with N F < 20, N B < 8 
in Ref. 0. The most important result from microscopic 
calculations 5 - is the prediction of instability regions, spe- 
cially when the number of 4 He atoms is small; neverthe- 
less, a core with five or more He atoms is able to bind 
any number Np of 3 He atoms. 

Our purpose is to analyze the ground state and the 
low-lying excited states of the 3 He atoms in a mixed he- 
lium droplet, by using the DMC method, and relying 
on the single-particle orbital orderings obtained previ- 
ously^ from the study of a single fermion. The resulting 
spacings resemble those of the rotational spectrum of a 
diatomic molecule, where the 4 He core plays the role of 
one atom and the 3 He is the other atom. At low energies 
each level is uniquely classified by its angular momentum, 
but vibrational-like excitations appear at higher energies 
for heavy enough clusters. The adopted filling order for 
large number of 3 He atoms is therefore Is lp Id . . . , 
with some attention being paid to a possible 2s intruder. 
We shall concentrate on clusters with eight and twenty 
bosons, varying the number of fermions. This moderately 
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large number of bosons ensures the existence of several 
bound excited levels. All calculations have been made 
using the HFD-(B) potential of Aziz and co-workers^ 
for the He-He interaction. 

The plan of the paper is as follows. In Section II we 
give some general details concerning the DMC calcula- 
tion. Section III collects single-fermion results needed 
later. Section IV deals with the case of two fermions: 
this system is simple enough to be analyzed in depth for 
a large set of quantum numbers, shedding light in partic- 
ular on the single-particle ordering. Section V is devoted 
to stable mixed clusters, for the specific cases of max- 
imum spin and maximum orbital angular momentum. 
The case of fully polarized clusters is also considered. In 
Section VI we exhibit the calculated one- and two-body 
distribution functions, both for normal and fully polar- 
ized mixed clusters. Finally, in Section VII we summarize 
our findings and draw some general conclusions. 

II. THE DIFFUSION MONTE CARLO METHOD 

The DMCSiSi method is based on an importance sam- 
pling wave function, with the double role of controlling 
the variance of the ground state energy and to incorpo- 
rate both the required statistics and the desired spin and 
angular momentum quantum numbers. In this work we 
use the same form employed previously^ The trial wave 
function is written as the product of four factors, 

*T = $BB ( f>BF ( S>FFD F (1) 

corresponding to a Jastrow form for the boson-boson 
(BB), boson- fermion (BF), and fermion-fermion (FF) 
parts, and the Slater determinantal part Dp for the 
fermions. Each Jastrow term is a translationally invari- 
ant and symmetric wave function, with the structure 

®MN=J[e^»^ (2) 

where indices M, N represent bosons (B) or fermions (F) , 
and indices i,j run over the corresponding atoms. The 
product runs over all different pairs i < j if M = N, and 
no restrictions apply when M ^ N. We have used a sim- 
ple but nevertheless physically complete representation 



The last factor Dp in the trial wave function is the 
determinantal part, which will be discussed later on for 
each specific case. The description of the trial wave func- 
tion is completed by the inclusion of the Feynman-Cohcn 
backflow 25 in the fcrmionic exponential tail as well as 
in the radial dependence of the Slater factor. Following 
Ref. 12(1 we have replaced each fermion coordinate by 
a transformed coordinate 

v i =v i + ^r 1 {r ij ){r l -Y j ) (4) 

For the backflow function rj{r) we choose the medium- 
range form used in Ref. |53 

fl(r) = A ( 5 ) 

keping the same value for the parameter A = 5 A 3 . 

Using the model described above for the importance 
sampling guiding function, a real-time DMC evolution 
has been carried out employing a 0(r 3 ) approximation 
to the Green function, with time slice r = 0.00025 K _1 , 
for a total of 200 blocks of 200 steps each. A block av- 
erage was used in order to diminish the unavoidable cor- 
relations of the DMC method, with the aim of obtaining 
a reasonable estimate of the variance. An initial pop- 
ulation of 1000 walkers lead, on the average, to a total 
number of forty millions samples. 

In spite of this block averaging, the resulting estimate 
of the variance turned out to be quite optimistic, i.e., 
much smaller than reasonably expected. As a conse- 
quence, we opted for a very costly but safe procedure, 
by carrying out ten independent calculations with the 
same time slice and the same number of samples, but 
with randomly selected initial set of walkers. The result- 
ing ten results lead to the true variance (or so we hope). 
The numbers quoted in what follows correspond to this 
prescription. 

Due to the presence of the Slater determinants, the im- 
portance sampling function is not definite positive. In our 
calculations, the random process has been constrained 
using the so-called fixed node approximation: any walker 
attempting to cross a nodal surface is neglected. As it 
has been shown^SiiSi the use of this approximation leads 
to an upper bound to the lowest energies. 



fMN{r) = —^y—-\ t pain , (3) 

depending on three parameters, b, p and v, in general dif- 
ferent for each of the three BB, BF and FF pairs. The 
short range coefficients b and v have been fixed, indepen- 
dently of the number of bosons or fermions, to the values 
v = 5.2 for any subset MN, bsB = 2.95 A, DBF — 2.90 A, 
and bFF = 2.85 A. The long range part rpMN entering 
Eq. ©, has the role of confining the system and fixing 
roughly its size. The three parameters pmn are adjusted 
for each droplet by minimizing its ground state energy. 



III. THE ONE FERMION SYSTEM 

The system made of a single 3 He atom plus a drop of 
4 He atoms has been recently investigated^ For the sake 
of completeness we summarize in this section the most 
relevant results for the analysis of systems with more 
fermions. 

As there is only one fermion, the Slater determinant in 
Eq. becomes a single-particle wave function 4> n e, m (r) 
of radial quantum number n and orbital angular momen- 
tum i with projection m. Translational invariance is en- 
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sured by referring the fermion coordinate to the bosonic 
center of mass. 

Different values of I select a specific angular momen- 
tum subspace in which the DMC procedure will drive the 
wave function so as to minimize the energy. Thus, the 
radial part of the fermionic wave function can be taken 
to have the form 



4>Um{r) = r l Yi m (f), 



(6) 



which ensures that there are no radial nodes (n = 1). 
These single-particle wave functions do not contain a ra- 
dial confining term, because it is already included in the 
Jastrow part of the importance sampling wave function. 
Notice that within the DMC procedure it is not possi- 
ble to obtain excited levels with n > 1, because this will 
require to impose a strict orthogonality on the physical 
n = ground state. Nevertheless, one may use an indi- 
rect procedure based in the moment method which will 
provide an upper bound to the energy of the radial ex- 
cited levels. 

Table H] displays the energy of these one-fermion states, 
as a function of the single- fermion quantum numbers. 
The quoted energies are slightly different from those 
of Ref. because of the improved statistics of the 
present calculation. The ground state energies of the 
pure bosonic droplets are also given, so as to define the 
separation energies 



E n t(N B ,N F = 1) - E{N B ,N F = 0), (7) 



IV. THE TWO FERMION SYSTEM 

Adding two 3 He atoms to a core of 4 He atoms results 
in a system deeply resembling the helium atom. To a 
large extent, the bosonic sub-cluster plays the role of the 
atomic nucleus, with the two 3 He atoms corresponding 
to the electrons. There are two families of levels, sin- 
glet (S — 0) and triplet (S = 1). Each of the states 
is characterized by the configuration, the orbital angular 
momentum L, the spin S and the parity P. 

Specific two-fermion states are constructed from the 
single-particle wave functions © by coupling the angu- 
lar momentum part to the required quantum numbers 
and by symmetrizing the singlet or antisymmetrizing the 
triplet radial wave functions. In general one does not ob- 
tain real wave functions, but the remedy is simple. For 
M = the result is real. Otherwise construct two cases 
for M and — M and either add or subtract them. The 
wave function thus constructed has no good L z , but still 
good total orbital angular momentum, and should have 
no effect on the computed energy. 

For the sake of completeness we list in Table [H] the 



TABLE II: Space part of importance sampling singlet(upper 
part) and triplet (lower part) wave functions for the two- 
fermion drop. 



L Config * (SINGLET) 



also quoted in the Table. Positive values refer to unbound 
levels. The energy obtained for the n = 2,1 = excited 
level is above the dissociation limit for Nb — 8 or very 
close to it for N B = 20. 



TABLE I: Energies of droplets with eight and twenty bosons 
and one fermion, for several values of the angular momentum. 
The row labelled Nf — is the system without fermions. 
The columns labelled e are the separation energies. Values in 
italics correspond to unbound levels. 



Config 


N B = 8 


N B = 20 




E(K) e(K) 


E(K) e (K) 


N F =0 

Is 

lp 

Id 

2s 


-5.14(1) 
-6.08(1) -0.94 
-5.60(1) -0.46 
-4.98(1) 0.16 
-5.10 0.06 


-33.76(2) 
-35.55(1) -1.79 
-35.15(2) -1.39 
-34.55(2) -0.89 
-33.80 -0.04 
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The single-particle wave functions of Eq. JBJ will be 
used later to construct model wave functions for systems 
with two or more fermions. Unfortunately the moment 
method does not provide a wave function for the radial 
excitation, and for the calculations of the following sec- 
tions we will use the simple form 



=o(r) 



(8) 



without any radial node. The lack of nodes may be a de- 
ficiency, but the truly important point is that the simple 
form chosen is linearly independent of the Is state. 
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specific forms thus obtained, and used as importance- 
sampling wave functions to drive the DMC stochastic 
procedure. 

TABLE III: Binding energies, in K, for several states with 
two fermions and Nb = 8 (upper table) and Nb = 20 (lower 
table), classified by the configuration and angular momentum 
quantum numbers. Values in bold face correspond to the 
physically interesting states. The dissociation limit is 6.08 K 
(Nb = 8) and 35.55 K (N B = 20). The statistical errors of 
the energies are between 0.01 and 0.02 K. The row labelled 
Eff. is the prediction of the non-interacting fermion model. 



N B =8 
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lsld 
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lp2s 
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7.09 
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7.02 


6.54 
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5.92 


5.44 


4.82 


6.02 




N B = 20 


L F S 


Is 2 


lslp 
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lsld 


lpld 


Id 2 
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0+0 


37.33 
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Eff 


37.34 


36.94 


36.54 


36.44 


36.04 


35.44 


35.59 





The energies obtained from these configurations are 
displayed in Table lilTI for N B = 8 and N B = 20. Some 
among these results lie above the dissociation limit (i.e., 
the energy of the lowest bound state with the same num- 
ber of bosons but with a single fermion) and do not cor- 
respond strictly to truly bound systems. 

It has to be understood that in the results quoted in 
these two tables, only a single number per row does have 
a physical sense: the DMC algorithm improves system- 
atically the importance sampling wave function, but be- 
cause of the use of the fixed-node approximation, that 
improvement only provides a variational upper bound to 
the energy for each of the subspaces with well defined L, 
S and P quantum numbers. To give an example: the 
state 1 5 l may have projections on the Is 2 , lp 2 , Id 2 and 



FIG. 1: The two- fermion spectrum for Nb = 8 and Nb = 20, 
classified accordingly to the configuration. Energies are in 
Kelvin. The horizontal line is the dissociation limit. Note 
that apart from the energy shift, the scales of the two plots 
are the same. 




Configuration 




Configuration 



other shell-model states, but the mixing will probably 
not be constructed along the DMC stochastic procedure. 
Thence, the DMC physically relevant results are those 
with larger value of the binding energy for each row. The 
other configurations with smaller energy are possible in- 
teracting configurations, and presumably an optimized 
linear combination within each row could provide a bet- 
ter binding energy. 

Note that among the states involving the 2s shell none 
of them has maximum binding except ls2s : 3 S because 
the configuration is unique. For Nb = 8 this level is not 
bound, but it is the last bound one for Nb — 20. For the 
other cases the DMC optimized energies prefer Is over 
2s orbits, but ls2s : 3 S survives because Is 2 : 3 S violates 
the exclusion principle. 

Concentrating on the normal shells (Is, lp, Id) it ap- 
pears that the binding energies depend basically on the 
configuration, and are almost independent of the cou- 
pling within the configuration. Figs. and [21 give an idea 
of this near independence, both for energies and radii. 
Within a central-field shell-model description, this fact 
indicates that the residual interaction between the two 
3 He atoms is very small, with energies close to the ones 
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FIG. 2: The root mean square radius (in A) of bosons and 
fermions for two-fermion systems, referred to the full center- 
of-mass of the drop, as a function of the configuration. The 
two plots correspond to Nb = 8 (upper figure) and Nb = 20 
(lower figure) . The lines displaying an almost horizontal lines 
represent the boson radii, and the varying lines represent the 
fermion radii. 




Configuration 




Configuration 



provided by the simple non-interacting picture 

E(ni£ 1 ,n 2 £ 2 ) = E Q + e nill + e n2 i 2 , 

where E is the energy of the boson core and e n £ the 
separation energies defined in Eq. Q and given in Table 
[I] The corresponding values appear in the last row of 
Table llTTl with the exception of the last column, the \p2s 
configuration which actually has evolved close to the lslp 
configuration. 

The rough picture of the two-fermion drops as basi- 
cally non-interacting fermions bound to a rigid 4 He clus- 
ter will be refined in Section IVII bv introducing an effec- 
tive monopole interaction. 



V. CLUSTERS WITH MORE THAN TWO 
FERMIONS: BINDING ENERGIES 

In this Section we consider selected states with a fixed 
number of bosons and an increasing number of fermions, 
up to 18 for normal clusters, i.e., with fermions with spin 
up and down, and up to 9 for fully polarized clusters, with 
only spin up fermions. 



A. The Slater determinant 

The fermionic factor Dp must have good angular mo- 
mentum and spin quantum numbers, and it must be 
properly antisymmetrized and translationally invariant. 
We take it to be a product of two Slater determinants, 
one for each spin orientation. The obvious way to build 
them up is with single-particle functions 4>n£m generated 
by a central field which dictates a natural filling order. In 
general, good total orbital angular momentum L and to- 
tal spin S demand a linear combination of determinants. 
If we insist on a single product, the construction is quite 
cumbersome for shells with high £ (see e.g. Ref. l28|) . 
though simple when dealing with s and p shells. 

The main problem in establishing a reasonable shell- 
model description of a system containing 3 He atoms, 
both in pure as well as in mixed drops, is the lack of phc- 
nomenological information about the central field. There 
are two familiar schemes common to other fermionic 
systems: the shell ordering Is, 2s, 2p, 3s, 3p . . . , charac- 
teristic of atoms, and the harmonic oscillator sequence 
Is, lp, (2s, Id), (2p, If) . . . , used in light and medium 
atomic nuclei. The principal quantum number follows 
different rules in both schemes. For the latter the paren- 
thesis indicate degenerate orbitals. 

In the previous Sections we have explained why the 
Is, lp, Id, 2s, . . . ordering was adopted for the calcula- 
tions. The information comes from the analysis of a 
single 3 He atom bound to a medium-size bosonic drop, 
which has been studied by density functional meth- 
oc jgI2ii4 and by microscopic methods based either in vari- 
ational wave functionsi^ii or in DMC techniques^ All 
these studies suggest that the fermion may be viewed as a 
particle bound to a potential well, centered in the bosonic 
drop, but which is appreciably different from zero (and 
attractive) only in a rather wide region near the surface 
of the drop, but goes to zero near the center of the drop 
as well as at long distances. This single-particle potential 
gives a special level ordering based on the orbital angu- 
lar momentum, Is, lp, Id, If ... , with almost degenerate 
single-particle energies. Moreover, the fermion has a very 
small probability penetrating the boson drop. The same 
scheme results from the study of one 3 He atom in liquid 
4 He, giving rise to the so-called Andreev states, as well 
as from the study of many 3 He atoms attached to a large 
core of 4 He atomsiSS Note, however, that intruder lev- 
els, like 2s, may appear for a sufficiently large number of 
fermions and a sufficiently small number of bosons. 

Among the possible wave functions related to a given 
configuration (the so-called terms in Atomic Physics) we 
have chosen the two simplest cases: a) maximum total 
spin on which maximum total orbital angular momen- 
tum is built; b) maximum orbital angular momentum on 
which maximum spin is built. The resulting wavefunc- 
tions are products of two determinants, one for each spin 
orientation. Other choices demand linear combination of 
Slater determinants. 

As the exponential tail in the trial wave function has 
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the role of roughly confining the system, we can construct 
the Slater determinants considering only the angular mo- 
mentum part of the single-particle functions, as well as 
the spin part. As mentioned above, we have used the har- 
monic polynomials (f>£ m (r) — r^l^ m (fi) as single-particle 
functions. The determinants so constructed are transla- 
tionally invariant wave functions, in the sense that they 
only depend on the 3Np — 3 relative coordinates — r 3 . 
This fact is of crucial importance, particularly when de- 
scribing systems with a small number of constituents. 

The way of constructing the required determinants is 
very simple. Take for instance the maximum spin case. 
Once the innermost shells have been filled, the remain- 
ing 3 He atoms occupy the </>«, </>«_i . . .spin up states, 
until the angular momentum states are exhausted; then 
the same procedure is followed to fill out the spin down 
states. The spin S of the resulting determinant has the 
maximum value allowed for the occupancy of the shell, 
and its orbital angular momentum is L = \L Z \. In general 
the value of the determinant is a complex number, not 
very adequate for the DMC algorithm. The solution is 
as in the two-fermion case: to use the sum or difference 
(whichever is non vanishing) of the determinants with 
L z = L and L z = —L. The importance sampling wave 
function has well defined S, Sz and L, but not L z . Nev- 
ertheless, due to the rotational symmetry of the Hamilto- 
nian, this has no influence on the energy values. However 
the need of computating determinants with complex ma- 
trix elements still remains, with the consequent slowing 
down of the numerical calculations. 

By using this procedure we have calculated states with 
up to Np — 18, corresponding to the complete filling of 
the three lowest shells. We have also considered states in 
which all spins are up, i.e., fully polarized fermions, with 
the maximum number Np — 9. The procedure to con- 
struct the determinants in this case is an obvious adap- 
tation of the one described above. 

In previous work£*£ we used a conventional Cartesian 
ordering, in particular {x 2 , y 2 , z 2 , xy, xz, yz}, whereby 
the single-particle orbitals are a mixture of 2s and Id 
wave functions. As a consequence the differences with 
the present calculations — associated to changes in the 
nodal structures — become significant when the d— shell 
starts to be filled. 



B. Binding energies of normal clusters 

Table IIVI presents the values of the binding energies 
corresponding to two situations, Nb = 8 and Nb = 20, 
for values of Np from to 18. All these values have 
been computed with the level ordering discussed above, 
and for two possible couplings: S max , where in each shell 
particles are aligned to maximum spin S, and then to 
to maximum orbital angular momentum L; L max , where 
particles are first aligned to maximum L, and then max- 
imum S. In Table llVl whenever there are two entries for 
a given cluster, the upper row corresponds to S max , and 



the lower one to L max . Maximum spin is quite uniformly 
favoured but the splitting of the two computed levels is 
always smaller than 0.3-0.4 K. 

TABLE IV: Binding energies (in K) of mixed clusters with 8 
and 20 bosons and up to 18 fermions. Whenever there are two 
entries for a given cluster, the upper row corresponds to the 
Smax coupling: particles aligned to maximum 5* and then to 
maximum L. The lower row is for the L max case: maximum 
L first and then maximum S. In the last column the results 
obtained in Ref. |g for Nb = 8 are also displayed. 
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The present results for Nf < 8 should coincide, and 
indeed they do, with previous calculations based in the 
cartesian ordering of the single-particle states^ Beyond 
Np = 8 the calculations of Ref. used an uncontrolled 
mixture of Id and 2s states, as explained above, thus cor- 
responding to different importance sampling functions. 
The last column of Tablc llVl displavs the results obtained 
in Ref. [U for Nb — 8 clusters. One should keep in mind 
that both results are based on the DMC method within 
the fixed node approximation, so that in both cases the 
obtained energies are actually upper bounds to the real 
ones. The present binding energies for Np > 8 are 
slightly higher that the previous results, the gain being 
of 0.16 K for Np = 9 and monotonically increasing up 
to 1.40 K for Np = 18, in the case of Nb = 8. This 
apparently modest increase (up to 6%) may be relevant 
for the boundaries of the stability chart of mixed drops. 
In any case, it gives support to the level ordering used 
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in the present work, as the associated importance sam- 
pling function provides a better variational bound than 
the previous ones. 

FIG. 3: Binding energies (in K) for Nb = 8 and Nb = 20, as a 
function of the number of fermions Nf- The long horizontal 
lines correspond to normal clusters, and the short lines to 
polarized clusters. Note that apart from the energy shift, the 
scales of the two figures are the same. 



have been plotted in Fig. [31 together with those for polar- 
ized clusters [to be discussed in Section IV Dl (Tabic IVT)] . 
The scale is such that values for the normal clusters are 
superimposed. 

Basically the energies grow linearly with the number 
of fermions but finer details emerge when looking at the 
fermion chemical potential 

tx F (N F )=E(N B ,N F -l)-E(N B ,N F ), (9) 
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FIG. 4: Separation energies, or ionization potentials, (in K) 
for normal S m ax clusters as a function of the number of 
fermions Nf- The lower curve correspond to 8 bosons, and 
the upper curve to 20 bosons. Dashed lines are just an eye 
guide. 
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which is plotted in Fig. 0] for two values of Nb and 
N F < 18. Observe the sudden drop of fi F after Np = 2 
and 8, corresponding to the closure of the Is and lp 
shells. The relative minimum appearing at Np — 14 is 
at first somewhat puzzling. Its origin will become clear 
in Section IVTl 



The binding energies of Table IIVI for normal clusters, 



C. The stability map revisited 

As it has been mentioned in the introduction, one of 
the most appealing properties of the mixed He clusters is 
the existence of instability islands, namely regions around 
selected values of (Nb,Nf) in which the system is not 
bound. These regions were discovered after Many-Body 
computations based on self-adjustable variational func- 
tions constructed with Jastrow factors supplemented by 
2p-2h and 3p-3h Configuration Interaction correlations.- 
Afterwards, the calculation was refined by means of the 
DMC method, 8 confir ming the previous findings. Given 
that both calculations provide really with only upper 
bounds to the energy, and having observed the improve- 
ment of the present DMC approach, based upon a differ- 
ent ansatz for the determinantal part of the importance 
sampling guiding wave function, we have revised the pre- 
vious calculation, just to check and eventually improve 
the limits of the instability regions. 

Indeed, in the previous subsection we have seen that 
the present level ordering leads to a noticeable energy 
gain in the ld-shell with respect to previous works. This 
fact suggests that some of the clusters previously qual- 
ified as metastable, i.e. systems with negative energy 
but less bound than clusters with a smaller number of 
fermions, could be in fact stable. Indeed, a new computa- 
tion near the beginning of the l<i-shell indicates two new 
bound systems, namely the clusters (Nb — 3,Np = 11) 
and (Nb — 4, Np — 9). Special attention has been paid 
to the cluster (Nb — 1,N F — 18), corresponding to full 
Is, lp and Id shells, but our finding is that this cluster 
is not bound. 

The stable clusters are displayed in Fig. \5\ which su- 
persedes the results previously obtained in Ref. 0- From 
the experimental point of view, we suggest the study of 
the regions Nb < 4 and Np < 5 to ascertain the stability 
limits. The measurements will require an improved mass 
resolution, at least 1 amu for clusters up to 25-30 amu. 
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FIG. 5: Stability map of mixed clusters. Solid squares 
represent truly bound states, and open squares represent 
metastable states. 
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D. Polarized mixed clusters 

We have also considered fully polarized mixed systems. 
The values obtained for their binding energies are shown 
in Table Ivl for Nb = 8 and 20. In both cases the number 
of fermions Np < 9 is limited by the complete filling 
of the Is, lp and Id orbits. The asterisk in Table El 
corresponds to situations in which the binding energy of 
the fully polarized system E^(Nb, Np) is smaller than 
the energy E(Nb, Nf — 1) of the normal cluster with one 
fermion less. In the other cases the results correspond to 
true bound states. 

The energies of polarized clusters are compared with 
those of normal clusters in Figure [3J Polarized mixed 
drops are always less bound than the unpolarized cluster 
with the same number of fermions. In other words, they 
are excited states. 

The energy differences between normal an polarized 
clusters take values around 0.25 K (N B — 8) and 0.15 K 
(Nb = 20) per fermion. It is worth mentioning that the- 
oretical calculations for liquid 3 He provide a difference of 
around 0.10 K per particle at the equilibrium density and 
—0.10 K at densities close to the solidification onei^i In 
other words, the preferred phase would be the polarized 
one at high densities. This anomalous behavior has been 
interpreted as a side effect related to improper nodal sur- 
faces for the unpolarized systems. The particle density 
of the fermionic phase in our case is much smaller than 
that of the fermionic liquid, and in consequence we can- 
not ascertain if our shell- model filling scheme will present 
such an anomalous crossing at higher densities. 

In the same figure one may appreciate a sudden change 
for the differences Ap = E DOlmEl i(Np) — E po i(Np) at 
Np = 4, which has a simple interpretation: it corre- 
sponds to the filling of the lp shell for the polarized case 
and the related jump in the ionization potential after 
adding a new fermion. 

As for the unbound clusters, all of them have very high 
spin, and their binding energy is larger than the energy of 
the polarized cluster with one fermion less. Therefore it 
is very likely that the set of polarized states corresponds 
to a stable branch above the dissociation limit, analogous 
to the so called displaced terms in Atomic Physics. 

The separation energies for the polarized clusters, dis- 
played in Fig. El follow basically the same pattern than 
the separation energies for the normal clusters. Again 



TABLE V: Binding energies (in K) for fully polarized clusters, 
for 8 and 20 bosons. The 3 He single-particle configuration 
is indicated in the second column. Results marked with an 
asterisk correspond to clusters with a binding smaller than 
the cluster with a fermion less according to Table HV1 



AT 


Config. 


T 

Li 


b 


AT Q 
iV_B = 8 


AT on 













5.14(1) 


33.76(1) 


l 


Is 1 





1/2 


6.08(1) 


35.55(1) 


2 


IsHp 1 


1 


1 


6.65(1) 


37.05(2) 


3 




1 


3/2 


7.30(1) 


38.61(2) 


4 


ls 1 ^ 3 





2 


8.06(1) 


40.22(1) 


5 


lsHp^d 1 


2 


5/2 


8.26(1) * 


41.37(2) 


6 


IsH/ld 2 


3 


3 


8.53(1) * 


42.57(2) 


7 


lsHp^d 3 


3 


7/2 


8.96(2) * 


43.83(3) * 


8 


ls x lp z ld A 


2 


4 


9.53(2) * 


45.19(3) * 


9 


lsHp 3 !^ 5 





9/2 


10.19(2) * 


46.63(2) * 



they recall the atomic ionization potentials with the sud- 
den drop once a given shell is closed. Note that for po- 
larized fermions the closure of shells occurs at Np = 1 
(Is), Np = 4 (lp) and Np = 9 (Id). 

FIG. 6: Separation energies for polarized clusters. 
2.00r 
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EFFECTIVE MONOPOLE INTERACTION 
ANALYSIS 



In this section we shall treat the results of our calcula- 
tions as data, and assume the validity of the shell-model 
scheme to find an effective one- plus two-body Hamil- 
tonian that could reproduce them. This two-body part 
will only refer to fermions as the boson cluster will be 
assumed as a fixed core that generates the single particle 
energies in Table [U leading to a one body potential 



U 



where n s is the number of particles in shell s. The sum 
is extended to the occupied shells. 

The two-body part is defined by matrix elements 

V r L s % = {rs:(LS)\V\tu:(LS)) 
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where \tu : (LS)) is a two-particle state in shells t and 
s, coupled to orbital angular momentum L and spin S, 
properly antisymmetrized and normalized. 

To simplify matters, the two-body part of the full 
Hamiltonian will be separated into monopolc (to) and 
multipole (M) contributions H = H m s + Hm- The 
monopole Hamiltonian H m s is defined by the property of 
giving the average energy of configurations at fixed num- 
ber of particles n r and spin S r . The closed shells and the 
one-particle and one-hole states built on them are con- 
figurations with a single state (we call this set cs ± 1). 
Hence, their energies are entirely given by H m s- For the 
other Slater determinants entering our calculations H m 
gives an average value that will be split by the multi- 
pole term Hm- We shall assume (and check) that the 
influence of Hm is small, and simply neglect it. 

The extraction of effective interaction averages goes 
back in time, and an important earlier reference is the 
work of FrenchiS The form we shall use here came later 33 
and was used to describe shell formation properties in nu- 



The counter-terms in the second line now ensure its van- 
ishing at the closed shell as well as at one-particle and 
one-hole states. As a consequence, their energies are fully 
given by the first line in Eq. H14fl . which we refer to as H m 
from now on. The advantage of this operation is that it 
decouples the determination of the V rs and b rs centroids, 
so that we can proceed with the former first, as they are 
the ones that give the global features. 

In principle, the six necessary centroids — V ss , V sp , V s d, 
V pp , V p d and Vdd — could be extracted from Table IIIII 
However, this parameter-free choice has large uncertain- 
ties and it is better to reserve it as a consistency check 
with the results of a more precise fit to the energies in Ta- 
ble [W] which we call Ei, i = 0, 18. It is very instructive 
to start doing the fit by hand, i.e., step by step. 

Upon filling, the closed shells become new "cores": 
E-2 = Ecs , E$ = Ec p ■ ■ ■■ The single particle energies 
arc taken from Table [Q Then 
E 2 = E Cs =E Q + 2e s + V ss . Extract V ss 



Ez — E, 



Cs 



2V S1 



E, 



Cs 



e p . Extract V sp 



clear physics^- A forthcoming review article^ contains To extract y we do not re i y on # 4 , because it is not 



comprehensive information on formal properties of H, 
and Hm- 

The average matrix elements are defined as 



purely given in terms of centroids, but on 
E 7 = E Cs + 5e 



lOVpp or 



V* 



Vj£.(2L + 1)(25 + 1)(1 + (-1) L+S 5 rs ) 



PV 



LS rsrs 



£ £S (2L+1)(2S+1)(1 + (-!)*+**„) 



Z L V r L s?s(2L + l)(l + (-l) L+s S rs ) 
£ £ (2L + 1)(1 + 



(10) 



where V rs is the full (scalar) average of two-body matrix 
elements, whereas Vf s are (vector) averages at fixed S. 
It is convenient to introduce the following combinations 



= - (3V r \ + V r ° s ), b rs = V r \ - V r ° s (11) 



(12) 



4 D r - 1 

where D r — 2(2l r + 1) is the maximum number of parti- 
cles in the shell. The standard result is then 

1 



H 



mS 



U - 



r<s 



n r (n s — S rs ) + 



3n r 



(13) 



where S r is the total spin operator corresponding to the 
particles of shell r, S r = ^2 ier cn/2, and the 3n r /4 sub- 
straction ensures the two-body nature of H ms by making 
the spin contribution vanish for single particle states. It 
has the drawback of producing a non-zero values at closed 
shells and single hole states. Therefore it is preferable to 
rewrite 



H mS =U + 



1 



V (1 + S rs ) 
r<s 



V rs n r (n s - S rs ) + 



3n r (D r -n r ) x 
A(D r - 1) 



(14) 



E 8 = E Cp = E Cs + Qe p + 15Vj 

The fit becomes overdetermined, signalling a problem 
with some basic assumption about the effective interac- 
tion, which we shall try to identify later. As of now let 
us settle for a compromise value of V pp . The next step is 
Eg = E Cp + e d + 2V sd + QVpd = E Cp + Ed, that deter- 
mines V s d + 3V p d- As the two matrix elements will only 
appear in this linear combination, the number of param- 
eters is reduced to five. Finally, for Vdd we have the same 
compromise problem we had for V pp . To find reasonable 
values for V pp and Vdd it was decided to do an overall fit of 
the five parameters. In principle, the idea does not seem 
very sound because three parameters are apparently well 
determined. As we shall see, this may not be the case, 
and the numerical fit will turn out to be sound. 

The hand-made fit involves cs ± 1 states that are com- 
mon to the Smax and L max cases in Table Hvl In doing 
the numerical fit, only the S m ax states were included. 
The results for the binding energies in the Nb = 8 clus- 
ters are given in Fig.0 The agreement is quite excellent, 
but for the fully polarized case, also shown, there are sig- 
nificant discrepancies that can be cured by introducing 
the full H m s through a single parameter b rs = b so that 
the contribution of the second line in Eq. (|14f) becomes 
b[S{S+l)/2-J2 r 3n r (D r -n r )/A(D r -l)}. The results 
of the fit are given in Table IVII The fitted and calcu- 
lated curves become nearly undistinguishable. As noted 
at the end of Section IV Bl the full energies are rather 
smooth patterns that tell us little about details. As a 
first approximation, Fig. [7] for the normal S ma x clusters 
is reasonably well represented by a straight line, which 
would be the analogue of the famous Bethe-Weizsacker 
"Liquid drop" formula for nuclei. The truly sensitive 
quantities are the separation energies (chemical poten- 
tials) in Eq. © and Fig. And, indeed, the true test 
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FIG. 7: Comparison of DMC binding energies with the 
monopole Hamiltonian with spin (continuous line) and with- 
out spin terms (dashed line). Circles correspond to the com- 
puted DMC values for normal Smax clusters. The lower line 
corresponds to the fully polarized case, and the upper group 
to the normal clusters. In both cases is Nb = 8. 
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TABLE VI: The fitted centroids in K 

N B = 8 N B = 20 

Value Error Value Error 



V ss 0.073 0.006 -0.019 0.007 

Vs P 0.079 0.002 0.080 0.003 

V pp 0.081 0.002 0.031 0.003 

(Kd + 3Vp d )/4 0.078 0.001 0.045 0.001 

V dd 0.069 0.001 0.045 0.002 

b 0.071 0.005 0.063 0.007 



of the monopole description comes in Fig. |H1 A smooth 
linear approximation to the binding energies would result 
in a constant. By introducing H m there is an enormous 
improvement in that the shell effects at closures are well 
reproduced (dashed line). However, the more detailed 
pattern between closures demands the S(S + 1) term in 
the full H m s'- the agreement with DMC becomes truly 
quantitative (full line). 

The numerical fit was made for the normal case we 
have called S max , but Table II VI contains another nor- 
mal mode, L max . As was noted, the hand-made fit is 
the same for both couplings, and it gives results that 
are almost as good as the numerical fit for S max , and 
very good ones for L max . But here the numerical fit also 
does a slightly better job, shown in Fig. El where H m s 
is seen to reproduce beautifully the staggering pattern 
between S = and 5 = 1/2 states (referred to as L max 
case in Table IIV|) . Here we are faced with some inter- 
esting physics: the numerical fit was chosen to find good 
compromise values for V pp and Vdd, but it does slightly 
better than a fit restricted to those over-determined pa- 
rameters. The hint is that the other three parameters 
are not as well determined as the "hand-made" fit sug- 
gests. The most likely reason is to be found in size ef- 
fects: as fermions are added, the overall radius evolves, 
and for a self-bound system it should go asymptotically 

1 /3 

as Np . [A strong indication in this sense will be found in 



FIG. 8: Comparison of calculated separation energies with the 
monopole Hamiltonian with spin (H m s, continuous line) and 
without spin terms (H m , dashed line), for .Smax states. Circles 
correspond to the computed DMC values. The lower group 
corresponds to Nb = 8, and the upper group to Nb = 20. 
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Fig. ^J] Therefore, the effective matrix elements should 
also evolve in a way our simplified H rn s ignores: The nu- 
merical fit then emerges as the sound and natural way to 
define a best compromise value, not only for V pp and Vdd 
but for all the monopole parameters. 

FIG. 9: Same as Fig. |5]but for clusters with L ma x- 
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Finally, Table lvTTl compares the centroids V® s and V^ s 
obtained with the fit to the 19 normal S max clusters with 
the values obtained directly from the DMC calculation 
of two-fermion states, again with a very good agreement 
(within the large errors of about 0.04 K associated to the 
latter), thus confirming the consistency of the effective 
interaction interpretation. 

It appears that very hard DMC calculations lead to 
results amenable to a very simple and cogent interpre- 
tation in terms of the monopole Hamiltonians H m s- 
The S(S + 1) contribution is particularly interesting: 
Fig. [7] suggests the idea that for large enough number 
of fermions the polarized clusters could become ground 
states. Though this is only a speculation, it may be also 
be taken as a strong invitation to push the study of mixed 
clusters much further. 
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TABLE VII: The two-fermion centroids as obtained from the 
effective interaction (columns labelled Eff ) compared with the 
DMC computed values. 
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pp 
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36.47 


36.50 
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36.19 


dd 


35.54 


35.54 


35.61 


35.70 



VII. THE SHAPE OF MIXED CLUSTERS 

A. Normal mixed clusters 

In this section we present several figures related to the 
shape of mixed clusters. In Fig. ^| there are the values 
of the root mean square radii for bosons and fermions 
referred to the center-or-mass of the cluster, for the se- 
lected cases Nb — 8 and 20. There are some fluctua- 
tions, probably related to the use of the mixed estimator 
method to compute these radii, and thus depending on 
the quality of the importance sampling wave function. 
Apart from these fluctuations, the most noticeable prop- 
erties which emerge from these plots are the almost con- 
stancy of the bosonic radii and the smooth growing of the 
fcrmionic radii. This manifests clearly the representation 
of the cluster as a quite rigid bosonic core with an halo 
of fermions. 



FIG. 10: The values of the root mean square radii (in A) for 
bosons (squares) and fermions (circles) referred to the center- 
of-mass of the cluster, as a function of the number of fermions 
Nf in the cluster. The number of bosons are fixed to Nb = 8 
and 20. 
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the one-body distributions of bosons and fermions with 
respect to the center of mass of the mixed cluster are 
displayed. These distributions are given by 



N, 



p M (r) = (*|££(r-[r,-R])|*), 



where label M stands for B (bosons) or F (fermions), 
Nm is the number of atoms of the given species, and R 
is the center-of-mass of the full drop. Given that these 
distributions may have an angular dependence for open 
shells, we have computed the spherical average of them. 
The distributions are normalized to the number of parti- 
cles of a given species, 



p M (r)dr = N M - 



We observe in Fig. ^2 that bosons are located in the 
same central region, being slightly compressed as the 
number of fermions increases. This shrinking is more 
important for the light Nb = 8 cluster, and is almost 
negligible when Nb — 20. With respect to the distribu- 
tions of fermions, they are clearly located at the surface 
of the bosonic subcluster, with a small penetration near 
the center of the drop in the case of Nb — 8, more im- 
portant for larger values of the number of fermions. In 
the case of Nb = 20 the dominating picture is that of a 
rigid core of bosons with a fermionic halo. 

A complementary information about the shape of clus- 
ters is provided by the two-body distributions, 



N M 



PAf(r,r') 



N m (Nm - 1) 



-W^r-r^r'-r,)!*), 



This picture is confirmed by the plots of Fig. 1111 where 



which is normalized to 1. 

Because of the finite size of the system under consid- 
eration, this distribution function depends on two coor- 
dinates, r and r', or, equivalently, on the distance of the 
center-of-mass of the pair ((r + r')/2) to the center-of- 
mass of the system and the relative distance (r — r') of 
the two particles, thus producing a function very difficult 
to plot. In order to get a more friendly quantity we have 
averaged the above two-body distribution with respect to 
its center-of-mass and for the remaining dependence we 
have computed the spherical average. The reduced pair 
distribution so obtained is now normalized to 1, 



pi 2 {r)dv = 1. 



B. Polarized mixed clusters 

There are noticeable similarities between the distribu- 
tion functions related to polarized systems and those cor- 
responding to normal clusters. The subcluster of bosons 
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FIG. 11: The density distributions (in A -3 ) of bosons and FIG. 12: The boson-boson and fermion-fermion distributions 
fermions with respect to the center-of-mass of the cluster, for (in A" 3 ), for the two selected cases N B = 8 and 20. 
the two selected Nb = 8 and 20 cases. 




Distance of atoms to cm. (A) 




Distance of atoms to cm. (A) 



is again hard, its radius being practically independent 
of the number of spin-aligned fermions Np. Differences 
with respect to the unpolarized cluster appear when com- 
paring the root mean square radii of fermions, as shown 
if Fig. I13l where the fermion halo in the polarized clusters 
is larger than in the unpolarized case. 

Something similar happens with the density distribu- 
tions of bosons or fermions with respect to the center 
of mass. The former remains basically unaltered when 
the number of fermions grows, and the latter follows the 
same pattern that in the case of unpolarized clusters. 

There is, however, a remarkable fact in connection with 
the two-body distributions and, specifically, the fermion- 
fermion distributions. These distributions are shown in 
Fig. El f° r the two selected cases N = 8 and 20. As can 
be seen there, the fermion-fermion distributions are very 
different from those obtained for the normal systems: the 
rise near 4 A is much less pronounced and the range is 
much larger, suggesting that the fermions are either less 
correlated or subjected to a long-range correlation. 

In order to appreciate the differences in the density 
distribution functions for normal and polarized clusters, 
we have plotted them in Figs.EDfor Np — 9 and Nb = 8. 
Apart from a clear change in the shape, the peak of the 



FIG. 13: The values of the root mean square radii (in A) 
referred to the center-of-mass, for bosons (filled squares) and 
fermions (filled circles for normal clusters and open circles for 
polarized clusters) as a function of the number of fermions 
Nf in the cluster. The number of bosons is fixed to Nb = 8 
(left panel) and to 20 (right panel). 




Number of fermions 



Number of fermions 



normal case is close to 4 A, while for the polarized cluster 
it is close to 7 A. 



VIII. SUMMARY AND CONCLUSIONS 

This work has been devoted to a detailed analysis of 
clusters made of 4 He and 3 He atoms. This new analysis 
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FIG. 14: The fermion-fermion distributions (in A -3 ), for 
Nb = 8 and unpolarized (top) and fully polarized (bottom) 
clusters. 
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^ 0.0002 - 




Atom-Atom distance (A) 



FIG. 15: Comparison of the fermion-fermion distributions 
(lower panel), for clusters with Nb = 8 and Nf = 9. The 
dashed line corresponds to the polarized cluster and continu- 
ous line to the normal cluster. 



H 0.0000 




5 10 15 20 
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was motivated by the improved knowledge of the order- 
ing of single-particle orbits, obtained after a systematic 
study of the spectrum of a single 3 He atom bound to a 
core of 4 He atoms^ which is expected to provide the 
optimal importance sampling trial function for the DMC 
calculation. The study has concentrated on clusters hav- 
ing a sufficiently large number of 4 He atoms so as to of- 
fer a simplified pattern: the bosonic constituents arrange 
themselves as a quite rigid core whereas the fermionic 
atoms are distributed in the surface of the bosonic sub- 



cluster, producing a halo. This arrangement, previously 
obtained by means of density functional methods, has 
been confirmed and pushed down to systems with a small 
number of constituents. 

One of the primary aims was to check the previously 
obtained stability map^ after the optimization of the im- 
portance sampling function, as well as the improvement 
of the fermionic nodal surfaces. No significant change 
occurred, and large instability islands are still predicted 
for a small (Nb < 3) number of 4 He atoms. 

The determination of correlation functions, partic- 
ularly the fermion-fermion distribution functions, as 
well as the analysis in terms of an effective interaction 
model suggests that the residual interactions between the 
fermions is very weak, of the order of the computational 
precision achieved (near 0.1 K). This fact is also reflected 
in the insensitivity of the energies to the (L, S) quantum 
numbers, the spectrum being essentially determined by 
the configuration. 

In addition to the normal fermionic phase, with a small 
value of the spin, we have also analyzed the possibility of 
having a ferromagnetic phase, with all spins aligned. The 
normal phase is energetically favored but the ferromag- 
netic one gives rise to bound states, even for a moder- 
ately large number of 3 He atoms. For example, a cluster 
with 20 bosons is able to bind up to six fermions, but 
beyond that number the system is above the dissociation 
limit. Nevertheless, because of the large values of the 
spin for the ferromagnetic phase, one may expect these 
unbound states to be long-lived, like in the case of polar- 
ized liquid 3 He, and thus to be experimentally detectable. 
Perhaps by circulating bosonic clusters through a cold at- 
mosphere of polarized 3 He atoms one could create these 
spin-aligned states sticking 3 He atoms one-by-one to the 
bosonic seed. 

The mixed systems may have a very rich excitation 
spectrum, because of the gap between the normal and the 
polarized phase. For example, for Nb = 20 and Np = 6 
there is a difference of 1 K between the normal state (con- 
figuration ls 2 lp 4 ) and the polarized state (configuration 
lsH^ld 2 ). The analysis of the intermediate filled con- 
figurations, such as ls^p^d 1 or ls 2 \p 2 ld 2 puts heavy 
demands on the DMC algorithm but it would be rela- 
tively simple in the density-functional method. Though 
it does not seem possible with the present experimental 
techniques to measure the spectrum, it is worth remem- 
bering that the excitation spectrum plays a relevant role 
in the production abundances of clusters / 3 -^ 3 ™! 

Finally one should stress the power of the effective 
monopole interaction analysis that points to the basic 
simplicity of the DMC results and invites to take on the 
challenge of unearthing the deep reasons of such simplic- 
ity. 
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